Correlation dynamics between electrons and ions in the fragmentation of D2 molecules 

by short laser pulses 



m 

O 

o 



Oh: 



X. M. Tong,* Z. X. Zhao and C. D. Lin 
J. R. Macdonald Laboratory, Physics Department, 
Kansas State University, Manhattan, KS 66506-2601 

We studied tlie recollision dynamics between the electrons and ions following the tunneling 
ionization of D2 molecules in an intense short pulse laser field. The returning electron coUisionally 
excites the ion to excited electronic states from there can dissociate or be further ionized 
by the laser field, resulting in D"*" + D or + D^, respectively. We modeled the fragmentation 
dynamics and calculated the resulting kinetic energy spectrum of to compare with recent ex- 
periments. Since the recollision time is locked to the tunneling ionization time which occurs only 
within fraction of an optical cycle, the peaks in the kinetic energy spectra provides a measure of 
the time when the recollision occurs. This collision dynamics forms the basis of the molecular clock 
where the clock can be read with attosecond precision, as first proposed by Corkum and coworkers. 
By analyzing each of the elementary processes leading to the fragmentation quantitatively, we iden- 
tified how the molecular clock is to be read from the measured kinetic energy spectra of and 
what laser parameters be used in order to measure the clock more accurately. 
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I. INTRODUCTION 

The fragmentation and ionization of D2 by intense 
optical laser fields has been an active area of theoret- 
ical and experimental studies during the past decades 
Q IE IE ^ IS ' In most of these experiments it was as- 
sumed that the D2 molecule is ionized in the early phase 
of the laser field producing ion which is subsequently 
ionized by the laser. Mechanisms for the ionization of 
ion include bond softening (SO) l a , charge resonance 
enhanced ionization (CREI) jl H llM El El El , in 
addition to direct ionization by the laser field. The dis- 
sociation and ionization of in the laser field result in 
D"*" -|- D or D"*" -I- D"*", with characteristic kinetic ener- 
gies refiecting the internuclear separation of the breakup 
of at the time when it is excited or ionized. Thus 
bond softening and CREI, which produce distinct peaks 
in the D"*" ion kinetic energy spectra have been observed 
experimentally and predicted theoretically. These peaks 
can be understood without reference to the ionization of 
D2 itself initially, i.e., the ionization of D2 and can be 
treated as two independent events. However, recent ex- 
periments H, IE 11,1s pointed out a new group of peaks in 
the D+ ion spectra at the higher energy (about 5 eV to 10 
eV per ion) which has now been attributed to the rescat- 
tering process 14, 15, 16, 17, 18, 19, 20, HI IH [H H . 
In the rescattering process, the electron which is released 
by tunneling ionization is driven back by the laser field 
to collide with the residual ion to ionize it or to ex- 
cite it. If the ion is excited, it can dissociate directly 
or be further ionized by the laser. In both cases, the 
D+ ion will have kinetic energy (the refiection principle) 
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characteristic of the internuclear separation where the 
ionization occurs. This paper is to examine all the ele- 
mentary processes that lead to the emission of D+ ions 
by the rescattering process following the initial tunnel- 
ing ionization of the D2 molecule. The calculation will 
be performed for D2 molecules only, but clearly the same 
model can be applied to H2 with minor modifications. 

In Sec. we first discuss all the elementary processes 
that lead to the dissociation or ionization of D2 molecules 
in a laser field. Starting with the tunneling ionization of 
D2, we address the following issues: (1) Calculation of the 
ionization rates of D2 from its equilibrium distance us- 
ing the molecular tunneling ionization (MO-ADK) theory 
[25|: (2) The classical trajectory of the ionized electron 
in the laser field and the Coulomb field of the ion, 
with initial longitudinal and transverse velocity distribu- 
tions following the description of the ADK theory [23| : 
(3) The free propagation and spreading of the nuclear 
wave packet after the tunneling ionization of D2 from 
its equilibrium distance; (4) Semi-empirical formulae for 
electron impact excitation cross sections of from the 
CTg ground state to the first few excited electronic states, 
in particular, the first (t„ and 7r„ states. These cross 
sections have to be evaluated at all values of internuclear 
separations and for different alignment angles of with 
respect to the laser polarization direction; (5) Evaluation 
of tunneling ionization rates of from the excited CTu 
and 7r„ states at each internuclear separation; (6) Follow 
the time evolution of dissociation and ionization dynam- 
ics to extract the kinetic energy spectra of the fragmen- 
tation products. While rates or cross sections for each 
of these elementary processes have been formulated, in 
the full calculation we only consider D2 initially aligned 
perpendicular to the direction of the laser polarization. 
In Sec. mil the resulting kinetic energy spectra of D+ 
are compared to the experiment of Niikura et al. Q 
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where D"*" ions were detected without knowing whether 
the other fragmentation product is a D or a D+. In con- 
tradiction to the conclusion of this work where the main 
peak in the D"'' kinetic energy spectrum was attributed 
to the dissociation of foUowing excitation to the cr„ 
electronic state, we conclude from our calculation that 
the main peak is due to the further ionization of the ex- 
cited by the laser. This has the consequence that 
the molecular clock we read is at a different time from 
the one read in Niikura et al. jj]. We further analyzed 
the contributions of the total kinetic energy spectra of 
D+ resulting from the different excited electronic states, 
from dissociation or ionization, and from rescattering at 
the first return, the third return, or higher returns. The 
simulated D"*" kinetic energy spectra from ionization also 
were compared to the recent experiment of Alnaser et al. 
|26j| . From our analysis, we conclude that the fragmenta- 
tion of D2 can be used as a molecular clock based on the 
rescattering dynamics. The clock can be read more ac- 
curately if the laser pulse is chosen at the lower intensity 
and with a shorter duration. With such a clock, the time 
duration can be read with an accuracy of fraction of a 
femtosecond without the attosecond laser pulses. Equiv- 
alently this means that the distance between the two 
nuclei can be read with an accuracy of fractions of an 
Angstrom. This can be achieved experimentally by com- 
paring kinetic energy spectra for experiments carried out 
at different mean laser wavelengths, or by comparing the 
kinetic energy spectra of D"*" and H"*" from the fragmen- 
tation of D2 and H2, respectively, in the same laser pulse. 
We finish this paper with a summary and conclusion in 
Sec.HVI 



II. THE THEORETICAL MODEL 
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FIG. 1: Schematic of the major physical processes leading 
to the formation of D"*" ion by the dissociation or ionization 
of D+. The first row depicts the oscillating electric field of 
the laser. The second row shows the spreading of the vibra- 
tional wave packet after the initial tunneling ionization. At 
to the D2 molecule is ionized. The tunneled electron returns 
to rescatter with the Dj" ion at ti where it excites to 
the excited (t„ or 7r„ electronic states. The excited can 
dissociate along these repulsive potential curves, or further 
ionized by the laser to produce D"*" ions by Coulomb explo- 
sion, to produce characteristic kinetic energy spectra of the 
fragments. For more detailed description, see Text. 



A. The elementary processes 

The schematic of the physical processes leading to the 
fragmentation of ion following the ionization of D2 
molecule in an intense laser pulse is depicted in Fig. ^ 
The D2 molecule is first ionized at tg near the peak of the 
laser pulse, releasing an electron into the oscillating laser 
field. At tl, the electron is driven back to the molecu- 
lar ion, to excite the other electron in the ion to one of 
the higher excited electronic states, or to ionize it. We 
will be dealing with peak laser intensities such that the 
returning electron does not have enough energy to ionize 
DJ. In the meanwhile, the nuclear wave packet propa- 
gates from its mean internuclear distance Ro=1.4 a.u. at 
to to Ri = 1.6 a.u at ti. Thus the electron impact excita- 
tion probabilities of the molecular ion by the returning 
electron have to be calculated for with a vibrational 
distribution ii). This distribution is indicated in 

the second row of Fig. ^ Once the is in the excited 
state, represented by the curve labelled (T„ and 7r„ in 
Fig. n the D2 can dissociate directly to D"*" -I- D, or it 
can be further ionized by the laser at t'^ when the elec- 



tric field of the laser returns to its peak value. If the 
is ionized at t ^ , then it will fragment by Coulomb explo- 
sion to produce D+ + D"^ ions. The total released kinetic 
energy for such a "two-step" process can be calculated. 

The rescattering does not have to occur only at the 
first return time ti. Due to the attractive field from the 
molecular ion, the released electron can return to collide 
with the molecular ion at later times, i.e., after more 
than one optical cycle, following the initial ionization. 
For example, the return can occur at t2 and t^, in the 
second optical cycle, or at t4 and ts, in the third optical 
cycle, and so on. At these later times excitation and 
ionization occur at larger internuclear separations, thus 
the kinetic energies of the fragmented ion are smaller. 
In general the returning probabilities become small after 
three optical cycles. Following the general convention we 
call t2 the second return and ta the third return, etc... 

An important feature of the elementary processes de- 
scribed above is that the rescattering times t^ and the 
subsequent tunneling ionization time t^ are relatively well 
locked to the clock tg of the initial ionization of D2. Since 
tunneling ionization occurs only near the peak of the 
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TABLE I: Relation between the returning time and the av- 
erage nuclear separation for and Dj. 



return 


time (fs) 


< R> (a.u.) 










ti 


1.9 


1.8 


1.6 


t3 


4.3 


2.5 


2.1 


ts 


7.0 


3.0 


2.6 


t7 


9.6 


3.2 


3.0 



laser field, to spans only a fraction of an optical cycle. 
Similarly, ti and t^ are also restricted to within sub-fs 
accuracy. These precise clocks in turn define precise in- 
ternuclear separations. For laser pulses with mean wave- 
length at 800 nm, the mean internuclear distances Ri for 
ti (1=1,3,5,7) are shown in Table|l]for H2 and D2. Note 
that at ty, the center of the vibrational wave packet for 
H;^ has already bounced back from the outer turning 
point, but not so for D^. A classical estimate shows that 
it takes 8.5 fs to reach its outer turning point. During 
these later times, the wave packet spreads significantly. 
The clock or the mean internuclear separation can be 
probed directly by the characteristic kinetic energy peaks 
of the fragmented D"*" ions. Changing the wavelength of 
the laser clearly will change the clocks and the mean 
internuclear separations. Replacing D2 by H2 will not 
change the clock but will change the mean internuclear 
distances. 

To read the clock from the measured kinetic energy 
of the fragmented ion, however, there are a number of 
factors that make the clock "fuzzy". First, the initial 
tunneling ionization occurs over an interval of about 0.3 
fs near the peak of the laser field. The initial vibrational 
wave packet, taken to be the ground vibrational wave- 
function of D2, according to the Frank-Condon princi- 
ple, has a width of 0.2 nm. This vibrational wave packet 
will broaden as it expands to larger internuclear sepa- 
ration. The electron impact excitation probabilities and 
the MO-ADK rates also depend on internuclear separa- 
tions. These factors would reduce the precision of the 
clock such that distinct peaks in the kinetic energy dis- 
tribution of the fragmented ions are not as clearly sepa- 
rated. We model the rescattering process to check how 
accurately the molecular clock can be read from the ki- 
netic energy spectra of the fragmented ions for different 
laser parameters. 

We now describe the models used for calculating the 
rates and probabilities for each elementary process. 



B. Tunneling ionization rates for molecules 

We first discuss how the ionization rates of D2 and 
in the laser fields are calculated. From the rescattering 
model above, we need the ionization rates for D2 from 
the ground state, and for from the excited CTu and 7r„ 



states over the whole range of R. The rates are needed 
for different alignment of the molecules as well. 

We calculated the tunneling ionization rates using the 
recently developed MO-ADK model '25'] . It was obtained 
by extending the widely used ADK _2X, theory for 
atoms in a laser field to molecules. In the MO-ADK 
theory the ionization rates are given in semi-analytical 
expressions. For a diatomic molecule in a parallel static 
electric field, the ionization rate for a valence electron is 
given by 



B'^{m) 



2I™I|to|! k22c/k- 



3 s 2Zc/K-|m|-l 



-2k^/3F 



with 



B(m)=^C,™(-ir 



'(2;-f l)(; + |m|)! 



2{l 



(1) 



(2) 



Atomic units are used unless 
Eq. QJ, K is related to the 



otherwise indicated. In 
ionization energy Ip by 
K = ^J'ilp, I is the orbital angular momentum of the va- 
lence electron, m is its projection along the internuclear 
axis, Zc is the effective charge seen by the valence elec- 
tron in the asymptotic region and F is the field strength. 
In Eq. |(2J), the parameters C;„i are determined from the 
valence electron wave function of the molecule in the 
asymptotic region. The laser peak power will be given 
in units of Iq = 10"'^'* W/cm^ and the mean wavelength 
is 800 nm. If the molecule is aligned at an angle with 
respect to the laser polarization direction, the ionization 
rate is given by 



m' 

where Wm' is given in Eq. except that 
B{m') = Y.^imDl.^^{0,e,Q) 



(3) 



x(-ir 



'(2^ + l)(;-f|m'|)! 
2{l-\m'\)\ ' 



(4) 



where the D-function expresses the rotation of the elec- 
tronic wave function from the direction of the molecular 
axis to the laser polarization direction. In the MO-ADK 
model, Eq. |^ reduces to the traditional ADK model for 
atoms if I is taken to be the orbital angular momentum 
quantum number of the valence electron. For diatomic 
molecules, the summation over / is a consequence of ex- 
panding the two-center electronic wave function in terms 
of single-center atomic orbitals. The coefficients Gim are 
functions of R and depend on the electronic states of the 
molecule. 

For D2 at the equilibrium internuclear separation, the 
parameters Cim have been calculated by Tong et al. [2^ . 
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FIG. 2: Binding energies of as a function of internuclear 
separation. 



Within the range of its ground vibrational wavefunction, 
it was found that the MO-ADK rates depend weakly 
on R. It was also found that the major component of 
B(m=0) in Eq. @ is /=0, for Da, thus the MO-ADK 
rates for D2 depend weakly on the alignment of the 
molecule. The accuracy of the MO-ADK rates for D2 
at the equilibrium distance had been checked previously 
and found to be in good agreement with the result from 
ab initio calculations p^ . 

We next consider the ionization of D^ in a laser field. 
Since the ionization rate depends sensitively on the ion- 
ization potential, in Fig.|21we show the electronic binding 
energies (R) at each R of the first four electronic states 
of D J . The negative of the electronic binding energy is 
the ionization potential. The total potential energy of 
each electronic state is Uj;(R)= Ei(R)+ 1/R. For peak 
laser intensity in the range of 0.5-5 Iq, estimate based on 
the simple ADK theory or the more complete MO-ADK 
theory shows that D^ in the ag state will not be ionized 
by the laser except for R greater than about 5 a.u., while 
for TTu and states the D^ will be readily ionized be- 
cause of the much smaller ionization potentials. Thus we 
need to calculate only the MO-ADK rates of D2 in the 
CTti state as a function of the internuclear distance R. 

In Fig. 131 we show the calculated MO-ADK tunneling 
ionization rates at F=0.06 a.u. for the (T„ electronic state 
of D^, for R in the range of 1-6 a.u. and for alignment 
angle 6*^0° and 90°. At 61=0°, the MO-ADK rates have 
been checked against the "exact" static tunneling ioniza- 
tion rates calculated usin g th e complex rotation method 
in the two-center system [S^. The MO-ADK rates tend 
to be somewhat higher, especially at small and large R 
region. For R greater than 6.0 the tT„ ionization energy 
is already very close to the ionization energy of atomic 
H, thus the ADK ionization rates of H(ls) are used for 
R > 6.0. In the actual calculation, the coefficients Cim 
are obtained for each R such that the MO-ADK rates 
can be readily calculated for any field strength, and any 



alignment angle of the molecule using Eq. 

C. The rescattering model 

Following the initial ionization of D2, a correlated elec- 
tron wave packet and a vibrational wave packet are cre- 
ated at to- The initial vibrational wave packet is taken to 
be the ground vibrational wavefunction of D2 , assuming 
that the ionization process is fast and the Frank-Condon 
principle is valid. Due to the heavy mass of the nuclei, 
the vibrational motion is not modified by the subsequent 
laser field. The time evolution of the vibrational wave 
packet is thus described by 

x{R,t) = ^ax.(i?)e-*^"*, (5) 

V 

a = y" Xg{R)Xv{R)dR. (6) 

Here {Xv{R)} and {et,} are the vibrational wavefunction 
and vibrational energy of D J in the ag ground electronic 
state, respectively, and Xg{R) is the the ground vibra- 
tional wave function of D2. 

The rescattering model for describing the motion of the 
ionized electron in the subsequent laser field is modeled 
similar to the method used by Yudin and Ivanov [2^ 
for He. The ionized electron is treated classically, under 
the combined force from the laser field and the residual 
Coulomb interaction from the D^ ion. For simplicity, the 
latter is approximated by an effective charge Zc=+1 at 
the midpoint of the internuclear axis. To calculate the 
trajectory of the ionized electron, we solve the equation of 
motion (Newton's second law), with the initial condition 
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FIG. 3: Static MO-ADK ionization rates for molecules 
aligned parallel (solid curve) and perpendicular (dashed 
curve) to the electric field direction. The filled squares repre- 
sent exact static ionization rates calculated using the complex 
rotation method I30II for parallel aligned molecules. The static 
field strength is 0.06 a.u.. 
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FIG. 4: Returning electron energy spectra for D2 in a pulse 
laser with peak intensity of 1.5 lo (lo^lO^'' W/cm^) and pulse 
length of 40 fs obtained from the simulation. 

that the ionized electron is at (x,y,z)=(0,0,zo) where zq 
is the tunnehng position from the combined potential of 
the Coulomb field and the static electric field. The initial 
velocity v is assumed to have a distribution from the 
ADK model, 

5(v) (X e--'«/^. (7) 

In this model, the tunneled electron is ejected isotrop- 
ically with a Gaussian distribution in velocity, i.e., we 
consider the ejected electron have initial velocity in both 
the transverse and the longitudinal directions. For each 
initial time to or phase (j)o that the ionized electron was 
born, the classical equation of motion was solved to ob- 
tain the trajectory. The distance of the electron from 
the center of ion is monitored for over seven opti- 
cal cycles for longer pulses or till the end of the laser 
pulse if the pulse is shorter. The distance of closest ap- 
proach of the electron from the ion and the time when 
this occurs for each trajectory are recorded. From these 
data, the impact parameter b and the collision energy T 
of the corresponding electron-ion impact (no laser field) 
excitation or ionization are obtained. 

Figure ^ shows the probability distribution of finding 
the returning electron with kinetic energy T measured in 
the asymptotic region for a laser with peak intensity at 
1.5 Iq. If the residual Coulomb interaction from the 
ion is neglected, the expected maximum returning energy 
will be 3.17 Up = 29 eV, where Up is the Pondermotive 
energy. The inclusion of Coulomb interaction increases 
this peak energy to about 35 eV. 

In Fig. ^ we show three groups of returning electrons. 
In the first group, the electron was born at to when the 
laser field has a positive phase (j)Q (i.e., beyond the peak 
field). It was driven outward and then back by the os- 
cillating laser field to recoUide with the ion within 
one optical cycle. This group is denoted by ti where 
the returning electron has peak current near 35 eV. The 
second group labelled t2 denotes an electron which does 



not collide with the ion at the first return, but at the 
second return about half a cycle later after the electron 
reverse its direction again. The kinetic energy for this 
group of electrons is smaller. The third group was de- 
noted by ta-l-t;^. For the ta, the recollision occurs at the 
third return. For the t 7 g roup, the electrons were born 
at a negative phase 0o |24j. i.e., before the laser reaches 
the peak field. These negative phase electrons do not 
recollide with the ion in the first optical cycle when the 
field change direction since they were accelerated by an 
increasing field right after birth. Due to the Coulomb fo- 
cusing by the ion they collide with the ion at the third re- 
turn. Without the Coulomb focusing the negative phase 
birth would not contribute to the rescattering process. 
In calculating the returning electron energy distribution 
shown in Fig. ^proper weights from the MO- ADK rates 
and the initial velocity distribution of the tunneling elec- 
tron have been accounted for. In Fig. ^we did not show 
the electron energy distributions from collisions occurred 
at returns after two optical cycles. The general trend is 
that at higher returns, the kinetic energy of the electron is 
smaller and the probability of rescattering is also smaller. 
In our calculations we have accounted for rescattering up 
to seven optical cycles for the long laser pulses. 

D. Electron impact excitation and ionization 
probabilities 

For each impact parameter b and kinetic energy T of 
the returning electron, we need to calculate the electron 
impact excitation and ionization cross sections of D^ at 
each internuclear separation R. Different from the He"'" 
case, there are few experimental or theoretical data avail- 
able for DJ. Thus we have to generate the cross sec- 
tions needed semi-empirically. For each total cross sec- 
tion cr{T) at kinetic energy T, we assume that the prob- 
ability for excitation or ionization at impact parameter 
b is given by 

P,n{b,T) = a{T) (8) 

ao = ^2//\E, (9) 

where T = /2 and A£' is the excitation or ionization 
energy. Here, the b-dependence is taken to be the Gaus- 
sian form. For the rescattering in He, Yudin and Ivanov 
[2^ have checked different forms of b-dcpendence and 
concluded that the results are rather insensitive to the 
precise functional form used. 

For electron impact ionization cross section, we employ 
the empirical formula 

a,{T,^E) = ^e'-'*^^^-'-'y^f{T/AE) (10) 
fix) = {Alnx + B{l--)-C—)-. (11) 

X XX 

where Ai? is the ionization energy. By fitting this for- 
mula to the accurate theoretical H(ls) ionization cross 
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FIG. 5: Ratio of the electron impact excitation cross section 
to 2po with respect to 2pi vs scaled excitation energy. Solid 
line is the fitted result from Eq. 112L open circles are from 
the calculation of Bray for H ^331 . and filled circles are from 
the He experimental measurement (3^. 



section ^ we obtained A = 0.7213, B=-0,302, C=0.225. 
The fitted formula, when applied to Hc+ , gives ionization 
cross sections in good agreement with the theoretical re- 
sults of Bray [3l| for He"*" as well. For at the equi- 
librium distance this formula also reproduces the recom- 
mended ionization cross section from NIST 32]. In this 
semi-empirical model, the molecular ion is treated as a 
point particle, thus the ionization cross section is inde- 
pendent of the alignment of the ion. 

For the excitation process, it is clear from Fig. |21 that 
(Tu and 7r„ states will be the dominant channels pop- 
ulated via electron impact excitation from the ground 



(Tg state since they have the lowest excitation energies. 
There are no theoretical or experimental data available 
for such cross sections as functions of internuclear separa- 
tions. Thus we will employ semi-empirical fitting proce- 
dure as well. We assume that the excitation cross section 
again can be fitted in the form of Eqs (|10|l and as in 
ionization, except that AE now is the excitation energy 
and the number 0.5 in Eq H10|l should be replaced by the 
excitation energy of the corresponding state in atomic 
hydrogen. From the tabulated H(ls)— > H(2p) excitation 
cross section in Bray [sj, we obtained A = 0.7638, B=- 
1,1759, C=-0.6706. The formula was further tested by 
comparing the predicted excitation cross section with the 
calculated one for e~-|- He"''(ls) e^+ He"'"(2p). From 
the total Is —>■ 2p excitation cross section, we can further 
distinguish excitation cross section to 2po or 2pi, with 
the direction of the incident electron beam as the quan- 
tization axis. The relative 2po and 2pi cross sections can 
be calculated theoretically or experimentally from polar- 
ization or correlation measurements. (Note: 2p_i cross 
section is identical to 2pi cross section by symmetry.) In 
Fig. O we show the relative cross sections of 2po to 2pi 
from the calculation of Bray [s^ for H, plotted against 
scaled energy (with respect to the excitation energy). On 



the same graph we display the same ratio for the excita- 
tion of He from Is^ to ls2p^P° from the experiment of 
Merabet et al. 34]. It appears that the both H and He 
data fall on the same curve when the collision energy is 
scaled with respect to the excitation energy. We fit the 
2po to 2pi cross section ratio by 



r{x) 



CTi 



8.2^1 + l.l/x^ 



0.44 



(12) 



where x=T/AE is the scaled kinetic energy. Since the 
ratio for He docs not differ much from the calculated 
ratio for H, this comparison convinces us to use the r{x) 
in Eq. (|12|l to describe the ratio for as well. The r{x) 
indicates that m=0 is the dominant magnetic component 
in the present interested energy regime. 

To relate the 2po or 2pi partial cross sections to the 
excitation cross sections of ct^ and tTu electronic states of 
D^, we need to know the alignment angle of the molecule. 
If the molecule is aligned along the laser field polariza- 
tion direction (which is also the direction of the electron 
beam), the 2po cross section is the excitation to ctu state 
and the 2pi (2p__i) cross section is for the excitation to 
the TTu state. If the molecule is aligned perpendicular to 
the laser polarization direction, then the role is reversed, 
i.e., 2pi (or 2p_i) corresponds to the cross section of the 
CT„ excitation, and 2po cross section to the 7r„ excitation. 
For any arbitrary alignment angle 9 of Dj, we assume 
the total excitation cross sections to (7„ and 7r„ are given 
by 



O-(crn) 

o-(7r„) 
ro 



(Tt {tq cos 9 + ri sin 
(TtItq sin^ 9 + ri cos 

(To + 2(Tl 

CTo _ r{x) 
ut t{x) -\- 2 
2(Ti _ 2 
gt r{x) + 2 



(13) 
(14) 
(15) 

(16) 
(17) 



The semi-emipirically fitted electron impact ionization or 
excitation cross section formulae discussed so far are for 
a free electron colliding with an atomic or molecular ion. 
For the rescattering process, the two electrons in D2 ini- 
tially are in the singlet state (S=0). Thus in principle, 
one should just use singlet excitation or ionization cross 
sections, instead of the spin-averaged cross sections. We 
obtain the singlet cross sections from the total cross sec- 
tion following the empirical formula derived in Yudin and 
Ivanov 23] [their Eqns. (8) and (9)]. 

These empirical formulae allow us to calculate electron 
impact excitation cross sections from Ug to Uu and to 
states at each internuclear separation and at each align- 
ment of the ion. In Fig. El we compare the electron 
impact excitation cross sections at the equilibrium dis- 
tance to (Tu and TTu states, for ions lying parallel and 
perpendicular to the incident electron direction which is 
also the direction of the laser polarization, respectively. 
When is aligned parallel to the laser polarization, im- 
pact excitation to Uu is the dominant channel. The Tr^ 
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to the 2s(Tg electronic state of D J . The empirical formula 
is chosen to be 



200 



50 100 
Electron Energy (eV) 
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FIG. 6: Electron impact excitation cross sections to au and 
TVu states of D J at the equilibrium distance, (a) The electron 
beam is parallel to the molecular axis; (b) the electron beam 
is perpendicular to the molecular axis. 



cross sections are smaller due to two factors: (1) the tTu 
state has higher excitation energy, see Fig. 12 (2) the 2po 
state has larger cross sections than 2pi for the electron 
energies considered, see Fig. |S1 The situation is different 
when ion is aligned perpendicular to the laser polar- 
ization direction. Fig. EJb) indicates that excitation to 
the TTu state is actually larger than that to the tJu state, 
at least in the 20-80 eV energy region. Note that in the 
experiments of Niikura et al. 0, Q the or were 
chosen to be perpendicular to the laser polarization di- 
rection. They assumed that electron impact excitation 
populates only the cr^ state, in disagreement with our 
analysis. 

The semi-empirical formulae presented above allow us 
to calculate electron impact excitation cross sections to 
(Ttj and TTu states averaged over the initially randomly 
distributed ions. We obtained the ratio of the cross 
section of cr,j with respect to 7r„, and compared the result 
with the ratio obtained by Peek js^ where the impact 
excitation cross sections for different internuclear sepa- 
rations were calculated using the Born approximation. 
The agreement is quite good, with the average cross sec- 
tion for (T„ about a factor of two larger than for 7r„. The 
absolute cross sections from Peek are larger since Born 
approximation was used. 

We also consider the small contribution from excitation 



<Je{T,AE) = 



1 



A£;2 

A 



f{T/AE) 
1 



1 + B/x X 



(18) 
(19) 



where the parameters A = 0.17, B=1.53 are obtained by 
fitting the formula to the Is 2s excitation cross sections 
of H. This cross section is assumed to be independent of 
the alignment of the molecular ion. 



E. Impact excitation probability 

With all the elementary cross sections available, we 
can now calculate the probability distribution of exciting 
D+ at a given internuclear separation R from the ground 
CTg state to a specific excited electronic state or to ioniza- 
tion states by the returning electron where the returning 
electron originates from the ionization of D2 molecule by 
the laser over a half optical cycle. The probability distri- 
bution is given by 



dPrn _ J J Pm{b,T)x^{R,tr)gMW{F cos (t,)dvd(t, 
dR ~ //3(v)W(Fcos(^)dvd(^ 



(20) 



The subscript m stands for the excited states (cr„, 7r„, ag) 
or ionization. Pm {b, T) is the impact excitation or ioniza- 
tion probability from Eq. 0. In this expression, W is the 
MO-ADK rate for ionizing D2 at the static field Fcos 0, 
where F is the peak field strength of the laser. For each 

0, the tunneled electron leaves the molecule with an ini- 
tial velocity v, with a distribution governed by Eq. ||7J), 

1. e., effects due to both the longitudinal and transverse 
velocity distributions are included. For each initial ve- 
locity and initial position of the tunneled electron, the 
return time tr at the distance of closest approach, the 
corresponding laser-free impact parameter b and kinetic 
energy T are calculated, and the excitation probability 
is also calculated. At each return time t,., the distribu- 
tion of the vibrational wave packet, x^{K,tr), is used to 
calculate the probability of finding at internuclear 
separation R. In this expression the MO-ADK rates and 
the impact excitation probabilities to ct„ and 7r„ states 
depend on the alignment of molecules. The other quanti- 
ties are isotropic. For D2 initially aligned perpendicular 
to the direction of the linear polarization of the laser, the 
impact excitation probabilities at different R's over half 
an optical cycle are shown in Fig.EI where the peak laser 
intensity is 1.5 Iq. Note that excitation probability to tTu 
is the largest, but to is also significant. On the other 
hand, excitation to 2sag excited state and direct ioniza- 
tion by the rescattering electron are not important. 

It is interesting to point out that the probability of ex- 
citation in Fig.|7|shows distinct sharp peaks as a function 
of R. To disentangle the source of these peaks, in Fig.|Sl 
we examine the contributions to the 7r„ excitation proba- 
bility according to whether the return time t^ falls within 
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one, two, three or four optical cycles after the tunneling 
electron is born. The excitation probabilities are larger 
for returns within one or two optical cycles. Within the 
first two optical cycles, the nuclear wave packets remain 
at small R with small spreading and the returning elec- 
tron has more kinetic energy (see Fig.^. For the higher 
returns the nuclear wave packet moves to larger R and 
spreads further and the smaller energies for the returning 
electron render the excitation probabilities smaller. 

We comment once again that with the inclusion of 
Coulomb attraction on the motion of the rescattering 
electron, the maximum returning electron energy is not 
given by 3.17Up=29 eV for the present peak intensity, 
but rather by 35 eV, as seen from Fig. 0] This has the 
effect of enhancing the excitation to the tt^ state as well. 

For peak laser intensity of 1.5 Iq the results in Fig. [3 
show that direct impact ionization of D J by the rescatter- 
ing electron is very small. The rescattering mostly popu- 
lates in the excited ttu and (t„ states. The dissociation 
of D J from an excited electronic state would release a to- 
tal kinetic energy given by U(Ro)-U(cxd), shared equally 
by D and D+, respectively. According to Fig.[7| excita- 
tion by the rescattering process peaks at characteristic 
internuclear separations related to characteristic rescat- 
tering time tr, thus measurement of the fragment ki- 
netic energies probes directly the recollision times. This 
forms the basis of molecular clocks in the experiments of 
Niikura et al. /J]- However, as shown in Tong et al. 
[36| and in Alnaser et al. 0|, the excited ions are 
still in the laser field and they can be further ionized by 
the lasers. Thus we need to calculate the kinetic energy 
spectra of D"*" resulting from Coulomb explosion after 
these excited ions are ionized by the laser. 




FIG. 7: Electron impact excitation and ionization probabil- 
ities of by the rescattering electron following tunneling 
ionizaiton of D2 by a short pulse laser with peak intensity of 
1.5 lo (Io=10" W/cm^) and pulse length of 40 fs. 




R (a.u.) 

FIG. 8: Electron impact excitation probabilities by the 
rescattering electron to the ttu state in the first four optical 
cycles after D2 molecules are ionized by a short pulse laser of 
peak intensity of 1.5 lo (Io=10^'' W/cm^) and pulse length of 
40 fs. 



Field ionization of excited ion 



In this subsection we consider the ionization of D2 
from the excited electronic states. We emphasize that 
we will consider peak laser intensity within 0.5-5 lo only 
where rescattering is important. In this intensity region, 
is readily ionized if it is in the n.^ excited state since 
its saturation intensity is only about 0.1 Iq because of 
small ionization energy. Thus we need only to calculate 
the ionization rate of from the cr„ state. If the ini- 
tial excitation to ctu occurs at R, the total accumulated 
probability for ionizing an electron by the laser field from 
the iJu state is 



P,(i?,(») 



^-fW{R')dt 

^-J^W(R')/v{R')dR'^ (-22) 



with 



-fiv^iR') = U{R)~U{R'). 



(22) 



where W{R') is the MO-ADK tunneling ionization rate 
described in subsection A, is the reduced mass of the 
two nuclei, and U{R) is the total potential energy of 
the <7u state. The cr„ state created at i?, followed by 
laser field ionization at R' will release a kinetic energy 
Ei{R') = U{R) - U{R') + 1/R'. Here we are more in- 
terested in the differential ionization probability which is 
given by 



dP.{R,R') ^ W{R') r^' wiR")/<R")dR" 



dR' 



v(R') 



(23) 



or in terms of differential probability per unit of kinetic 
energy 

dP,{R,R') dP,{R,R')dR' 



dE 



dR' 



dR' 
1 



dU 



\ dU{R') l 
I dR' I 



(24) 
(25) 
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FIG. 9: Kinetic energy spectra of ion resulting from 
laser ionization of D J in the excited state, for ion reaching 
the excited state initially at R=1.6 a.u. (solid curve) and 
R=2.1 a.u (dashed curve), respectively. Laser parameters: 
peak intensity 1.5 xlO^'' W/cm^, pulse length 40 fs. 



Fig. 1^ shows the expected ionization spectra from the (t„ 
state if it is initially created at two different R's (R=1.6 
and 2.1 a.u.), chosen to be the peak positions of the vi- 
brational wave packet at the first and the third returns. 
Clearly, the early return releases more energy (higher en- 
ergy peak). The spectra show many sharp peaks since 
ionization occurs only when the laser field is near its peak 
intensity at every half cycle. 

To obtain the total ionization spectra, we need to add 
up contributions from initial ionization at all values of R, 
i.e., 



dPic 



dE 



dPrndP,{R, R') 



dR 



dE 



dR 



(26) 



This integration is important primarily only for ioniza- 
tion from the excited (t„ state. For other excited elec- 
tronic states, due to the high ionization rate, ionization is 
complete within one cycle or less and we can set R = R', 
and the differential ionization spectra for these excited 
electronic states are given by 



dP, 



dPm dR 



dE 



dR dU' ^^"^^ 

The total ionization spectra are obtained by adding up 
contributions from all the excited electronic states, and 
from the initial ionization by the rescattering electron 
(very negligible). 

For the dissociation process, the energy spectra are 
obtained from 



dPr 



dis 



dE 



dPm dR 
~dR"dU' 



(28) 



The total dissociation spectra are obtained by adding up 
contributions from all the excited electronic states. In 
reality, the dissociation comes from the au excited state 
only. In all other excited electronic states the ions 
are immediately ionized by the laser in one optical cycle. 




10 15 
Energy (eV) 

FIG. 10: D"*" yield at several laser intensities for a 35 fs pulse 
length. The experimental data are from Ref. 4| for 1.5 lo 
where Io = 10^'* W/crn^. The peak values from experiment and 
from theory for 1.5 lo are normalized to each other. For peak 
intensities of 1.0 lo and 0.8 lo, the yields have been multiplied 
by 1.4 and 3.0 respectively to have the same peak ion yield 
height. 



III. RESULTS AND DISCUSSION 

The kinetic energy spectra of D+ ions can be deter- 
mined without any coincidence, as in the experiments of 
Niikura et al. |^ Q , or by detecting the two D+ ions in 
coincidence, as in the experiments of Staudte et al. 
or in Alnaser et al. |26j . we will present our simulation 
results for both types of experiments. 



A. Non-coincident D^ kinetic energy spectra 

In the experiments of Niikura et al. j^, 3; the ki- 
netic energy of D"*" ion was measured in the direction 
perpendicular to the direction of laser polarization. The 
measured D"*" signals come from ionization and from dis- 
sociation. Thus, 



Signal oc 2 



dP 

dE 



dE 



(29) 



In Fig. 1101 the experimental D+ kinetic energy spectra 
from Niikura et al. are shown. The energy scale is 
the total breakup energy, or twice the energy of the D"*" 
ion. The experiment was performed for a pulse of 40 fs 
and peak intensity of 1.5 Iq. We have shown simulations 
with the same laser parameters but with three peak laser 
intensities, at 1.5, 1.0 and 0.8 Iq. First we normalize the 
peak height at 12 eV between theory and experimental 
data at 1.5 Iq. Since the peak positions do not vary with 
laser intensity, we can normalize the calculated spectra 
at the two other intensities as well, with a multiplicative 
factor of 1.4 and 3.0 for the 1.0 and 0.8 Iq, respectively. 
If one compares the experimental spectra with the the- 
oretical one calculated at the same 1.5 Iq, clearly the 
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high energy peak near 16 eV from the theory is too high, 
while the theoretical spectra between 5 and 10 eV are 
somewhat too low. However it appears that the discrep- 
ancy can be reconciled if one takes into account of the 
volume effect in that the experimental spectra have to 
be integrated over a volume where the intensities are less 
than the peak value. The energy resolution and the finite 
acceptance angles can all contribute to the smoother ex- 
perimental spectra. One of course should also take this 
"better agreement" with caution in view that the peak 
intensity of the laser is often not known precisely. 

One of the major goals of the simulation is to unravel 
the origin of the structure in the kinetic energy spectra 
which in turn would provide insight of the working of the 
molecular clock. For this purpose, we show in Fig. 1111 
the calculated kinetic energy spectra, but separate the 
contributions from dissociation and ionization, and from 
rescattering occurred after one or two optical cycles, or 
equivalently, from the first (ti) or the third returns(t3), 
at two laser intensities, 1.5 and 0.8 Iq. At the higher 
intensity, in this figure we notice: (1) ionization is much 
stronger than dissociation; (2) the peak from the third 
return (2nd cycle) is higher than from the first return; 
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FIG. 11: Decomposition of D"*" ion yields into contributions 
from dissociation and ionization, and for rescattering occur- 
ring within the first and the second optical cycle after the 
initial tunneling ionization. The peak laser intensities are (a) 
I = 1.5 lo and (b) I = 0.8 lo, where Io=10" W/cm^ and pulse 
length is 40 fs. 



(3) The width of the peak from the first return is broader 
than the peak from the third return. The broadening is 
a consequence of the factor dR/dU in Eqs. (jSEl and ^ 
which is approximately given by R^. Another interesting 
observation is that the peak position of the dissociation 
spectra from the first return almost coincides with the 
peak position in the ionization spectra from the third 
return. This shift is due to the binding energy of the 
excited electronic states. 

In Niikura et al.'s experiment 4] the peak at 12 eV 
was attributed to originate from the dissociation of 
via the ct„ curve at the first return. In other words, this 
peak reads the clock at ti. According to our simulation, 
the peak comes from ionization following rescattering at 
the third return, and this peak should read the clock at 

t3. 

Contributions to the D+ signal from dissociation do be- 
come more important at lower laser intensity, as shown 
in Fig. Illf b'). Even at this intensity, the peak at 12 eV 
still comes mostly from the ionization following rescat- 
tering at t3 instead of dissociation following rescattering 
at ti. Furthermore the third return peak is higher than 
the first return peak for either dissociation or ionization. 
We remark that the spectra in Fig. ^| were calculated 
including contributions up to four or five optical cycles 
after the initial tunneling ionization and convergence of 
the calculation was checked. 



B. D coincident kinetic energy spectra 

The ion kinetic energy distributions in laser-D2 in- 
teractions have been determined in coincidence measure- 
ments where the two ions were detected simultane- 
ously by Staudte et al. and more recently by Alnaser 
et al. |2a|- In the latter experiment, the branching ratios 
of ionization with respect to dissociation had been mea- 
sured as well, for peak laser intensities of 1-5 Iq. Their 
data for peak intensity of 2.8 Iq are shown in Figure [T^ 
The experiment used a 35 fs pulse with mean wavelength 
of 800 nm. The spectra are from Coulomb explosion 
of ions at 60-80° with respect to the direction of the linear 
polarization of the laser field. In the figure we show the 
result of our theoretical simulation for laser intensity of 
2.0 Iq. We found best overall agreement with the experi- 
mental data at this intensity without considering volume 
effect and the fact that the theoretical calculation was 
carried out for molecules aligned perpendicular to the 
laser polarization while the experiments measured ions 
coming out of 60-80° with respect to the laser polariza- 
tion. The simulated spectra near the kinetic energy peak 
region of 7-12 eV agree quite well with the data, but the 
peak near 17 eV is more pronounced in the simulation. 

A direct comparison of simulated kinetic energy spec- 
tra with experimental data is complicated in general not 
only by the volume effect, the angular resolution of the 
D+ product, but also the difficulty of knowing the peak 
laser intensity precisely. In Fig. 1131 we show the yield 
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FIG. 12: Comparison of D+ ion spectra resulting from the 
double ionization of D2 molecules in a laser field. The exper- 
iment data are from Ref. (2^ for peak laser intensity of 2.8 
lo and the theoretical simulation is for laser peak intensity of 
2I0, where Io=10^* W/cm^ and the pulse length is 35 fs. 
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FIG. 13: Simulated D"*" ion yield from the double ionization of 
D2 by the rescattering process at several peak laser intensities 
in units of Io = 10^* W/cm^. The pulse length is 35 fs. 



for making two D"*" ions vs the total kinetic energy for 
peak laser intensity from 1.0-3.0 Iq. The calculations 
were done for 35 fs pulse and mean wavelength of 800 nm 
and with molecules aligned perpendicular to the laser po- 
larization direction. Clearly the yield increases rapidly 
with laser intensity. We further note that the peak po- 
sitions in the spectra do change with laser intensity. In 
particular, the main peak shifts to lower kinetic energy 
at higher laser intensity. To understand the reason of 
this shift, in Fig.^jwe separate the kinetic energy peaks 
into contributions from the ct„ and from the tt^ curves, 
and for rescattering occurring after one, two and three 
optical cycles following tunneling ionization. Recall that 
we consider D"^ from ionization only here. At 2.0 Iq, we 
note the larger contribution comes mostly from ioniza- 
tion of in the Tr^ state, although contribution from 



cr„ is not negligible. From Fig. I14f b) one can clearly 
identify the two peaks approximately at 10 and 12 eV 
in Fig. 1121 can be attributed to ionization from CTu and 
7r„, respectively, for rescattering collision from the third 
return. At 3.0 Iq (Fig. I14r a') ), due to the larger con- 
tribution from the CTu excited state, the peak positions 
in the kinetic energy spectra are shifted to lower values. 
Thus the sum kinetic energy spectra at the two higher 
intensities look different from those at lower intensities, 
as seen in Fig.^| Fig.^Ja) also shows that contribution 
from the third cycle becomes relatively more important 
at higher intensity. At higher intensity, the rescattered 
electron has larger kinetic energy. Thus it takes more 
time for the Coulomb attraction to bring the electron to 
come near the ion for the rescattering to occur. 

Figure ^] also illustrates how the working condition 
for using rescattering model to measure the precise time 
in a molecular clock can be limited. The kinetic energy 
spectra from each excited electronic state of D J has rela- 
tively well specified and distinct peak positions from the 
first, 3rd and 5th returns. Such peak positions imme- 
diately give information about the molecular clock since 
each peak position does not depend on the laser inten- 
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FIG. 14: D"*" yield from laser ionization via the (t„ (solid 
line) and tvu (dashed line) excited states with laser intensity 
(a) I = 3.0 lo and (b) I = 2.0 lo (Io=10" W/cm^). Each 
yield is further decomposed into contributions for rescattering 
occurring after one, two and three optical cycles. The laser 
pulse length is 35 fs. 
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sity. However, when ionization from (t„ channel also con- 
tributes then the combined sum would shift the peak 
positions as the laser intensity is changed, as shown in 
Fig. El Thus to read the molecular clock accurately, 
one has to choose laser intensity where only one of the 
excited electronic state contributes mostly to the ion- 
ization signal. Failure to do so would compromise the 
accuracy of the clock. Since the relative contributions 
of the ionization signals from cr„ and tTu are expected 
to change with laser intensity and with the alignment of 
the molecules, this also helps explain why the valleys in 
the experimental spectra are usually less sharply peaked 
than the ones simulated from the theory at a given peak 
laser intensity. 

C. Laser- H2 interactions and wavelength 
dependence 

Clearly the present method can be used to predict the 
kinetic energy spectra if H2 is used as the target. The 
only difference in H2 is that it has smaller reduced mass 
such that the wave packet propagates faster, and thus 
kinetic energy spectra will be shifted to lower energies. 
If the wavelength of the laser is increased, the period is 
longer and thus the kinetic energy spectra also will shift 
to lower energies. We have applied the present theoretical 
model to study the comparison of kinetic energy spectra 
taken for H2 and D2 simultaneously and also the 
variation of the kinetic energy spectra when wavelength 
was varied as in the experiment of Niikura et al. Q , see 
Tong et al. 

IV. SUMMARY AND CONCLUSION 

In this paper we have provided a comprehensive study 
of the elementary processes of the rescattering mecha- 
nism leading to the fragmentation of following the 
initial tunneling ionization of a D2 molecule in a short 
intense laser pulse. Ionization rates of from the ex- 



cited electronic states and impact excitation and ioniza- 
tion cross sections by the returning electron have been 
obtained based on the MO-ADK theory and from semi- 
empirical formulation, respectively. Following the initial 
idea of Corkum and coworkers, we showed that the ki- 
netic energy spectra of D+ in the higher energy region 
(5 to 10 eV per D+ ion) can be used as a molecular 
clock which can be read with subfemtosecond accuracy. 
Through our detailed simulation, we concluded that the 
dominant peak in the kinetic energy spectrum is due 
to the further ionization of the excited following im- 
pact excitation by the returning electron, and this exci- 
tation occurs not at the first return but mostly at the 
third return. We have compared our simulation results 
with the recent experiments of Niikura et al. and of 
Alnaser et al. with general good agreement. Further ex- 
perimental studies in terms of dependence on laser wave- 
length, pulse duration and alignment angles may provide 
more critical test on the present theoretical model. From 
the theoretical viewpoint, despite of the semi-empirical 
nature of the present modeling, we do not expect any 
meaningful pure ah initio quantum calculations viable in 
the foreseeable future. The present model has the further 
advantage that the mechanism for producing each indi- 
vidual peaks in the kinetic energy spectra can be iden- 
tified and the effect of laser parameters can be readily 
tested. On the other hand, the semi-empirical nature of 
the modeling can claim its reliability only after it has 
been exposed to more stringent tests from the experi- 
ment. 
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